arXiv:1505.04362vl [quant-ph] 17 May 2015 


The energy level structure of a variety of one-dimensional 
confining potentials and the effects of a local singular perturbation 

M.L. Glasser* 1,2 and L.M. NietcT 

'Dept. de Ffsica Teorica, Atomica y Optica, Univ. de Valladolid, 47011 Valladolid, Spain 
department of Physics, Clarkson University, Potsdam, NY 13699, USA 

May 19, 2015 


Abstract 

Motivated by current interest in quantum confinement potentials, especially with respect to the 
Stark spectroscopy of new types of quantum wells, we examine several novel one-dimensional singular 
oscillators. A Green function method is applied, the construction of the necessary resolvents is reviewed 
and several new ones are introduced. In addition, previous work on the singular harmonic oscillator 
model, introduced by Avakian et al. is reproduced to verify the method and results. A novel features is 
the determination of the spectra of asymmetric hybrid linear and quadratic potentials. As in previous 
work, the singular perturbations are modeled by delta functions. 


1 Introduction 

The bound state spectra in confining potentials, especially of linear and quadratic form, has been a concern 
in quark confinement, for example, for some time [T]. The effect on these spectra due to local changes in 
the potential at specific points is also of concern and led Avakian et al. [2] (see also mm to introduce the 
singular harmonic oscillator model V(x) = ax 2 + b5(x) which, along with generalizations, has accumulated 
an extensive literature. Of particular interest is the extensive body of fundamental work on these systems 
by S. Fassari et al. 001318], and that summarized in the book of Albeverio et al. [9] (see also uni) 
to which we refer the reader for further references. More recently, these models have arisen in the study 
of semi-conductor quantum wells where, since the 1980’s, it has been possible to engineer these entities 
with increasingly sophisticated properties. Quantum wells having parabolic confinement in applied electric 
fields are fundamental for a growing optical device industry. For a survey of the basic experimental and 
theoretical work in this area see, e.g. m- In the pioneering work these systems could be treated successfully 
as singularly decorated square wells and a Green function matching technique was developed providing an 
accurate account of the electro-luminescence of Ga-As based quantum wells m- In this note we present 
a first attempt to extend this work to the quantum-confined Stark effect in low symmetry and anisotropic 
structures now being explored. For details and references see m , e.g. 

Our aim in this note is to reproduce the basic energy level calculations for the harmonic case by the 
expeditious Green function method derived in [14, and then to go on to examine a number of new systems, 
including the important linear oscillator V(x) = a\x\ + bS(x — q) which, surprisingly, has not been studied 
in detail. For completeness, we begin by deriving the basic Green functions. 
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where 0(t) represents the Heaviside step function, is derived in many elementary texts, e.g. [T5| [16] and 
plays an important role in problems concerning electrons in magnetic fields, for example. However, its 
time-Fourier transform, the equally useful energy dependent Green function (Schrodinger resolvent) 


G ho (x, a/; = e iE+t ' h K ho {x, x'; t)dt 


(1.3) 


is less well-known, although it seems to have first appeared on p. 74 of the classical book of Titchmarsh 
El. In the case of the three-dimensional oscillator ( |1.3| ) was evaluated explicitly by Krebtukov and Macek 
[18 . In the next section we present an independent derivation of Gho(x , x '; E ) in one dimension in terms of 
the Parabolic Cylinder Function D v (x) (see Chap. 8 of USD, by obtaining an integral representation for the 
Hermite polynomials resolvent series. We study two applications of this result: the case of adding a linear 
potential to the harmonic oscillator and the case of the hybrid asymetric oscillator proposed in El; in both 
cases the spectrum is analyzed in detail. In section 3 we obtain the Green function for the linear confining 
potential V(x) = a 3 \x\, in terms of Airy functions. The applications to an asymmetric linear potential and 
to a couple of different combinations of oscillator and linear potentials are studied. In Section 4 we revisit 
the (Fdecorated harmonic well problem [2], and we analyze the same type of “decoration” for the linear 
confining potential previously mentioned. The paper ends with a summary and discussion. 


2 Calculation of the harmonic oscillator Green function 

In order to determine the Green function of a system uniquely, the boundary conditions must be “built 
in” to the solution of the Green function equation. In the case of (1.3) this is achieved by adjusting the 
integration contour. For the three dimensional oscillator, this was achieved [18j by assigning the energy an 
infinitesimal imaginary part. The subsequent integration was then carried out by taking advantage of the 
rotational symmetry of the system. In the one dimensional case this is not an option and we must proceed 
differently. 


From the familiar eigenstates ^°{x) and eigenergies E^° of the harmonic oscillator (1.1), as for any 
Sturm-Liouville system El, we form the Green function (resolvent) as a series, where the boundary 
conditions are automatically incorporated into the eigenfunctions, 
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As it is well known, this Green function satisfies the general equatior0 
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for the oscillator potential V(x) = muo 2 x 2 /2. Sometimes another dimensionless version of the Green 
function is used, let us call it G(x,x'] e = E/(faE)), related to G(x,x f ] E ) as 
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In this particular case, in order to simplify the calculations, instead of working with the function ( |2.l| ) 
let us set up the equivalent dimensionless system obtained by introducing the dimensionless quantities 
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which is a meromorphic function of s having non-positive integer poles. Thus if we obtain its value 
in any singularity-free region of the complex plane, analytic continuation can be invoked to determine 
its value elsewhere. Let us, therefore, begin by assuming that Re[ 5 ] > 2 and use the representation 
Jq 1 u n+s ~ 1 du = (n + s) -1 to get 
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where Mehler’s formula (see Sec. 10.13(22) of [19]) has been used to sum the series. Next, by means of the 
change of variables t = u 2 /(I — u 2 ) equation \2.1\ becomes the Laplace representation 
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After a simple calculation we obtain 

S(z,w,s) = e G 2 + w2 y 2 r(s)D_ s (V2 Z> ) D_ S (-V2 Z< ), (2.9) 

where z < (z > ) represent the smaller (larger) of z,w. The expression ( |2.9| is analytic throughout the s- 
plane, except at the poles of the gamma function. This expression for S does not seem to appear in the 
literature on Hermite functions. 


We can now replace in (|2.9|) z, re, s in terms of x, x', E and take into account Q2.5| )-(2.6) to obtain 

D e-l (pxy) D e _l , (2.10) 
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where now x< (x>) represent the smaller (larger) of x,x f , and the physical parameters /i, e are given in 
(2.3). The poles of ( |2.10 ), and therefore the bound states, are precisely the poles — n of the Gamma 
function, giving back the result 0- 

In Figure [l] we show plots of Gho(x, x f \ E) as a function of x and for a couple of values of the 
energy. One can appreciate a lot of symmetry in these plots due to the fact that the eigenfunctions of the 
oscillator entering the Green function ( |2.1| ) are real, and either even or odd, implying that G/ lo (x,x / ; E) = 
G/ lo (x / , x ; E) and Gho(x, x r \ E) = Gho {~ x ', E), as is clearly seen in the figures. 
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Figure 1: Green function ( 2. 10| ) 
and on the right Gho(x , x \ E — 


in units such that u — h — m — 1 for two energy values: on the left Gho{x , x \ E — 2) 
4.3). 


2.1 Harmonic oscillator plus uniform electric field 

Let us analyze first the case of a uniform electric field, given by the linear potential a 3 x, added to the 
original oscillator potential. By completing the square in the Schrodinger equation, it is easily seen that 
the system is again a simple harmonic oscillator, but with shifted coordinate x a = x + ip and shifted energy 
E a = huj(e — {fHf/2) 2 )\ 
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The Green function for this system is then simply 
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where all the physical parameters are given in ( |2.12| ). 

A plot of the Green function G a (x : x';E) is shown in Figure [2] for E = 2 and four different values of 
the electric field. The symmetry that appeared in Gho(%, x '\ E) is still present in G a (x, x'; E) because they 
are related through (2.13); the symmetry center is no longer the origin but the point (—<£>,—<£>) on the 
(x, x')-plane. The poles of (2.13), and therefore the bound states, are 


£ »°= r “(" + )-(?)’) 


n = 0,1,2 ,... 
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From these expressions it is quite obvious that in the limit a, <p —> 0 we recover the bound states energies 
and the Green function of the simple harmonic oscillator. 


2.2 Asymmetric harmonic oscillator potential 

As a second application, consider the asymmetric oscillator potential 
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Figure 2: Green function (2.13) in units such that cd = /i = ra=lfora fixed value of the energy (E = 2) and four 
different values of the electric held: from left to right, and from top to bottom, a 3 = 0,0.5, 1.3, and 2. 


and let Gj(x,x';E) denote the Green function G/ lo (x,x / ; E) given in (2.10) evaluated for uj = c Oj. Then 
the asymmetric oscillator states, which do not coincide with (n + 1/2 j = 1, 2, are given by [T4] 
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Note that, from equation (2.17), e(A) = e(l/A)/A, so we need only examine 0 < A < 1. For A = 1 (the 
symmetric harmonic oscillator = 0 J 2 ) there are clearly no new eigenvalues, apart from the already known 
(1.1). A plot of the eigenvalues e = E/(huj 1 ) as a function of A = uj\/uj 2 is given in Figure [3| where the 
harmonic oscillator eigenvalues are also indicated by the horizontal red lines. Note that in this plot it is 
quite clear that when A = 1, that is, for the symmetric oscillator, the correct eigenvalues are recovered. 

For A 0, that is, 002 oc and fixed, the potential (2.15) physically corresponds to an infinite 


barrier for x > 0 and a semioscillator of frequency uj\ for x < 0; in this case the spectrum is well known to 
be (2 n + 1 + 1 , something that is also clearly seen in the left part of Figure pj note that the states of 

the simple harmonic oscillator (A = 1) evolve as A —> 0 in such a way that 


3 7 

lime 0 (A) = -, lim e x (A) = 
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lim e n (A) = 2n + 1 + n = 0,1,. 
A—m 2 


(2.18) 


and then the spectrum e n (A = 1) = (n+|) is rescaled in a continuous way to become e n (A = 0) = (2n+l+^) 
in the limit A —0 or u ^2 oo. 

The case A —» oo can be imagined as UJ 2 0 and uj\ fixed, and corresponds to a semi-oscillator of 
frequency u)\ for x < 0 and V = 0 for x > 0, a situation without bound states, a fact that is suggested in 
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Figure 3: Bound states of the asymmetric oscillator e n = E n /(huji) vs A = uj\/uj 2 (solid curves), for small values 
of A (on the left) and bigger values of A (on the right). The harmonic oscillator energy eigenvalues are indicated 
with dot-dashed horizontal red lines, and the harmonic oscillator case with a vertical dashed orange line. 


the right part of Figure [3] by the fact that 

lim e n (A) = 0, n = 0,1 ,... (2.19) 

A—>-oo 

We should mention that, although in the present application we have concentrated in obtaining a 
closed transcendental expression for the bound states, it is also possible to obtain the Green function for 
the asymmetric potential. Nevertheless, the task is quite cumbersome and we reserve it for the future. 


3 Calculation of the Green function for V(x) 
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As a second example, let us consider now the Schrodinger equation for a linear potential of the form 
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A related problem was considered in Ref. eh- It is well known that equation ( |3.1| ) can be solved in terms 
of Airy functions. The general solution is 
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As we require ^(± 00 ) = 0, C 2 = 0. Now consider the differential equation defining the Green function 
G a (x,x'] E) for equation similar to ( |2.2| ), or G^(x,x';p) = —h 2 G a (x,x';E)/(2m), similar to ( |2.3| ): 


[dl + C 2 e ~ C 3 x'\ g) = 5(x - x'). 


(3-4) 


Clearly G^(x, x g) must be continuous for x — x’ and G^(x, x'\q) = G^(x', x; g). Hence, the solution must 
have the form 

G Ai((x — g) Ai(— (x' — p), x > x', 
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But, there is also the jump condition in the first derivative. In order to generate the ^-function on the 
RHS of (3.4), the partial derivative d x G^(x,x'] g) must be discontinuous: 
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In the present case we must have 


(3.6) 


CC [Ai((x' - 0) Ai'(-C^' - q) + Ai (~(x' - g) Ai'((x' - £>)] = 1 . (3.7) 


Since the quantity in (|3.7|) is essentially the Wronskian, it is independent of x', 


so we can set x' = 0 to get 


G a (x,x';E) = - 


f 2m\ 2/3 1 Ai(Ca;> - g) Ai(-Ca;< - g) 
V h 2 ) 2 a Ai(— g) Ai'(— g) 


(3.8) 


where, as in the previous section, x< ( x > ) represent the smaller (larger) of x, x', and the physical parameters 
g and ( are given in ( |3.3| ). The poles of ( |3.8| ) are the bound states of this problem: basically the alternating 
zeros of Ai \z) (even states) and Ai(z) (odd states). Formulas and list of numerical values of these zeros 
can be found in many references, for example [22 

In Figure [4] we show plots of G a (x , x '; E ) as a function of x and x', for a couple of values of the energy. 
As in the harmonic oscillator, a lot of symmetry is also present in these plots, a fact which is due to the 
symmetry of the potential, implying that G a (x,x f ; E) = G a (x',x; E) = G a (—x',—x;E), as it is clearly 
seen in the figures. 



Figure 4: Green function of equation (3.8) for h 2 = 2m and a = 1, that is, g = E and £ 
on the left G a =i(x,x'; E = 2.3) and on the right G a =i(x, x'; E = 3.2). 


1, for two energy values: 


3.1 Asymmetric linear potential 


As an application of the result we have obtained in (3.8), let us consider the asymmetric linear potential 


V(x) 


—afx, x < 0, 
a^x 2 , x > 0, 


(3.9) 


and let G 0) (x, xE) denote the Green function G 0 (x. x'; E) given in ( |3.8| evaluated for a = ctj. Then the 
asymmetric linear states, are given by the general result (|2.16|) (taken from jf4)): 


G ai (0,0;E) + G a2 (0,0- 1 E) = 0. 

In terms of g = (2m / H 2 ) 1 / 3 E/af and /3 = a\/a 2 , this can be written as 

_ Ai(—g)/Ai'(— q) = 

Ai(-gp 2 )/Ai'(-gp 2 ) P ' 


(3.10) 


(3.11) 
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A plot of this result is shown in Figure [5] There, it can be seen, for example, how the energy of the state 
n varies as ol^ changes, for aq constant: as aq 0 (/3 oo) the right branch of the potential goes to 
zero, and the bound states disappear; as ct 2 00 (/3 —> 0) the right branch of the potential becomes an 
infinite barrier and the discrete bound states that “survive” are only those coming from the solutions of 
Ai(— g) = 0 (those coming from Ai'(— g) = 0, which are also poles of (3.8), are no longer solutions of the 
limit case). 



Figure 5: Bound states of the asymmetric absolute value potential (3.9) given implicitly by (3.11): the physical 
parameter giving the energy, g — (2m/h 2 ) 1 ^ 3 E /of, is represented as a function of the frequency quotient /3 = olx/ol^- 


3.2 Half-oscillator-half-linear potential 

As an application of the results of the previous sections, let us consider an interesting example, the com¬ 
posite potential consisting of a harmonic oscillator for x < 0 and a constant force for x > 0: 


V(x) = - muj 2 x 2 6(—x) + a 6 x 0(x). 


(3.12) 


In this case, the equation providing the bound states ( |2.16| ) turns out to be 


G/io(0,0; E) + G a (0,0; E) = 0, 


(3.13) 


where G/ lo (0, 0; E 1 ) is given in ( |2.10| ) and G a (0,0;£’) in 


^ = (2mhuj 3 ) 1 ^ 6 /a = (2/(/up)) 1 
form 


Using the two parameters e = E/Tilo and 
where /i, p were already introduced in (2.12), equation ( 3.13| ) adopts the 


Ai'(-e 2 e) 
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r (i 
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r (!-f; 

In Figure [6] a plot of the first bound states of the composite potential ( |3.12 ) as a function of £ is shown. 
If a 0 (that is, £ —»• 00 ) the potential becomes zero for x > 0 and then the bound states disappear 
completely. If a —» 00 (that is, £ -> 0) an infinite barrier emerges for x > 0, and only the odd states of the 
harmonic potential survive. 

For the special value £ = \/2, a plot of the left hand side of equation (3.14) as a function of e is given in 


Figure [7J In this last plot the zeros of the function are clearly visible, and can be computed numerically. 
The calculation yields the results shown on Table [l] for the first ten bound states, which are graphically 
represented on the RHS of Figure [7| 























Figure 6: The first bound states of the asymmetric composite potential (3.12) from equation (3.14) as a function 
of £ = (2 mhu 3 ) 1 / 6 /a, for h 2 = 2m; the particular value £ =s y/2 is stressed with a vertical dashed orange line. 



Figure 7: For the particular value £ = >/2, a plot of the LHS of equation ( 3.14| (with ti 2 = 2m) as a function of 
e, where the zeros are clearly visible, is given on the left. On the right a plot of the numerical values of the zeros, 
given in Table [l] 


3.3 Harmonic oscillator and symmetric linear potential 

Let us consider now the following potential, which is clearly related to the previous examples: 

V(x) = ^ muj 2 x 2 + a 3 \x\. 


(3.15) 


Using the parameters e, /i, p and a defined in (2.12), the wave function solution of the corresponding 
Schrodinger equation, which is continuous at x = 0 and bounded at Too, can be expressed in terms of 
parabolic cylinder functions as 


1 p(x) = 


= C D a _ 1/2 (/J,x + fitp), x>0, 
foix) = C D <7 _ 1/2 (-iix + nip), x<0. 


(3.16) 


The results of [211 show us how to build the Green function for a problem like the previous one, gluing 
the solutions in two different regions (a, 0) and (0,5). In particular, as in the present case, if we consider 
a —>> — oc and b —» +oc, the Green function, as a solution of an equation similar to (2.3) or (3.4), is obtained 
taking the appropriate limits in equation (12) of [21], and we get 


G(x,x':E) = 


V>l(z>) V>2(^<) 

w[Mx'),Mx')Y 


G(x, x’\E) = — 


2m ^i(x>) V> 2 (ff<) 
h 2 W[tP 1 (x , ),Mx , )Y 


(3.17) 
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eo 

0.50501 


^5 

3.41789 

<T 

1.27615 


ee 

3.86844 

^2 

1.88901 


e 7 

4.29867 

^3 

2.43392 


es 

4.71332 

e 4 

2.94119 


eg 

5.11461 


Table 1: First ten bound states energies for the composite potential (3.12), evaluated numerically from 
(3.14) for the particular value £ = y/2. 


where the denominator in (3.17) is the Wronskian of the functions 'ipi(x') and ^(V). Using the two 
functions in (3.16), and evaluating the Wronskian at x' = 0, we get for the present problem 


_ 2 m + W) + w) 

nh 2 D ff _ V2 M [wAt-i/2(w) - 2 Ar+1/2 (w>)] ’ 


The poles of this Green function, that is, the zeros of the denominator, are given by 


A7-1/2 (m^) [ wDa - 1/2 (w) - 2-0^+172 (w)] = 0. 


(3.19) 




Figure 8: The bound states for the potential (3.15) obtained from equation ( |3.19| : on the left in terms of the 
parameters a and pup given in (2.12); on the right the dimensionless energy e = E/huo as a function of /up. The 
even states (orange curves) correspond to the solutions of {/J>(p) = 2D cr+1 / 2 (/tc^), and the odd states 

(blue curves) are the solutions of D a _ 1 / 2 (pp) = 0. 


It is possible to analyze the limit a,tp 0 (which implies a -A E/huo ) of this expression, which 


corresponds to the simple harmonic oscillator Green function (2.10) 


E) = \ j 


—r(-~ E 

Attijiuj V 2 hu) 
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i (/'•'•>) Dj L _i(-/ix < ) = Gho{x,x'\E). 


(3.20) 


4 Additional Dirac 5 interaction potential 

The appearance of the Dirac delta function and its derivative as potentials in problems of the sort we 
consider has been examined in depth by Albeverio et al. [9], whose work indicates that including them is 
compatible with the Green function techniques we have been using. In 1987 Avakian et al. |2j introduced 
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the singular oscillator model V(x) = \uj 2 x 2 + 0,8 {x — q), q = 0, and discussed its bound states. We shall 
generalize their results by letting q ^ 0: 


V a (x,q) = aS(x — q), a £ M. (4.1) 

This generalization will be studied for the unperturbed potential V u (x) to be both the harmonic oscillator 
and the absolute value potential considered in the previous section. 

We take into account the fact that the eigenfunction of eigenvalue E , corresponding to the potential 
V{x) = V u (x) + V a (x,q), must satisfy the Lippman-Schwinger equatiorj^] 

/‘OO 

if)(x) = - dx' G u (x,x'\E)Va{x',q)4>(x'), (4.2) 


with G u {x, x'\ E) the “unperturbed Green function” given either by (2.10) or by (3.8). 


4.1 Harmonic oscillator with additional 5 interaction 


Let us consider first the presence of a ^-perturbation, like in (4.1), on the harmonic oscillator potential. 
Then, the Lippman-Schwinger equation ( |4.2| ) is 

(4.3) 


ip(x) = -a dx' G ho (x, x'; E) 5(x' - q ) ip(x') = -a G ho (x, q; E) tp(q). 


Now, we set x = q and find that the eigenvalues of the energy are given by Gho{q , q]E) = —1/a, that is 

1 


tD € i {p) D e _i (-p) + 


rq-e) 


= 0 , 


(4.4) 


with e = E/hu> and p = pq. where p is given in (2.12) and r = . Note that (4.4) generalizes 

equation (12) of j2_. Indeed, for q = p = 0 our result completely agrees with that of Avakian et al. J2]. 



z 


Figure 9: First bound states of (4.4) as a function of the intensity of the Dirac delta perturbation r oc a, for two 
different positions of the singularity: p — 1/2 (blue line) and p — 2 (green line). The dotted red lines correspond to 
the energies of the unperturbed harmonic oscillator states, obtained also for r — 0 (orange dashed line). 


The result we have obtained is completely symmetric in g, therefore, as in [2], and without loss of 
generality, we will take in the sequel g,p > 0. In Figures [9| and [To| we show plots of the bound states 
of the oscillator modified by the presence of the Dirac delta perturbation V a (x,q). Note that if a,r < 0 
(attractive Dirac delta), the resulting bound state energies are lower than the unperturbed ones; on the 
contrary, if a, r > 0 (repulsive Dirac delta), the new energies are higher than the unperturbed ones. 

2 In operator form we have G U (H U — E) = I for general unperturbed Hamiltonian H u and Green function G u - If we 
consider an eigenstate |^) of the perturbed Hamiltonian H u + V a with eige nvalu e E, then (H u + V a — E)\ip) = 0. In addition, 
G U (H U — E)\ijj) = \ip) = —GuVal'ip)- In the ^-representation we have then ( |4.2| >. 
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Figure 10: Energy of the first bound states of (4.4) as a function of the singularity coordinate p = gq, for several 
intensities rocaof the Dirac delta perturbation. On the left we have three cases of “repulsive” delta-perturbations:: 
r = 0.5, (blue line), r — 1 (orange line) and t — 1.2 (green line); on the right, three cases of attractive delta- 
perturbations: r = —0.5, (blue line), r — — 1 (orange line) and r = —1.2 (green line). In both figures the dotted 
red lines correspond to the energies of the unperturbed harmonic oscillator states, obtained also for r — 0. Observe 
that as p grows, the perturbation on the lower states decreases. For either attractive or repulsive perturbation, at 
p — 0 the energy of the odd states do not change, but it does for the even states. 


4.2 Absolute value potential with additional S interaction 

Let us analyze now the effect of the presence of a S perturbation V a (x,q) = aS(x — q), a E M, on the 
eigenstates of the absolute value potential. This problem was considered in the past [23], but only for 
q — 0. In the present case case, the Lippman-Schwinger equation (|4.2|) is 


/ oo 

dx' G a (x , x '; E ) S(x' — q) !p(x') = —a G a (x , q\ E ) 'ip(q), (4.5) 

-oo 

where G a (x,x'; E) is given in (3.8). Now, we set x = q and find that the eigenvalues of the energy are 
given by G a (q,q; E) = —1/a, where G a (x,x'; E) in given in ( |3.8| ), that is, 

^Ai(Cg- q) Ai(— (q - g) = Ai(-p) Ai'(-p), (4.6) 

where g and ( are defined in ( |3.3| ), and 

2/3 


T] = 


a / 2m\ 


2 a \ h 2 J 


The complete Green function for this model can be obtained from the Dyson equation 

/ oo 

G a (x, y; E) V a (y, q) G(y,x’;E ) dy , 

-OO 

and turns out to be 


G(x, xE) = G a (x , x '; E) + « 


GJ^g^E) Ggis^x^E) 

1 -aG a (q,q;E) 


(4.7) 


(4.8) 


where G a (x,x'-,E) is given in ( |3.8| . 

In a similar way transcendental equations can be obtained giving the variation of the energy levels 
of the harmonic oscillator plus linear potential (2.11) and the harmonic oscillator plus symmetric linear 
potential (3.15) as a consequence of the addition of a Dirac delta perturbation as (4.1). 
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5 Summary and discussion 

In this report we have applied the Green function technique to study a number of simple, but apparently 
new, quantum mechanical problems dealing with confining quadratic and linear potentials in one spatial 
dimension. For completeness we have derived the Green functions involved, though at least one of them has 
been available in the literature for many years HSJ. Specifically, we have examined the harmonic oscillator in 
a constant force field, the asymmetric linear and quadratic oscillators, two hybrids thereof and two cases of 
decoration by Dirac delta potentials. By means of a composite Green function matching technique, we have 
streamlined the formulation of the transcendental equations determining the bound-state energy levels and 
their dependence on various parameters. We have treated the decoration problem by means of the Lipmann- 
Schwinger bound-state formula which leads to the necessary eigenvalue equations more expeditiously than 
by the wave function matching analysis used in prior studies mi of the decorated harmonic oscillator, 
confirming these pioneering calculations. We hope in the future to return to the construction of the exact 
Green functions for the composite oscillators. This work also serves as an exploratory study of models and 
techniques which we hope to apply to a new class of inhomogeneous quantum wells. Another natural step 
is to study stronger singular perturbations such as the Dirac delta derivatives considered in [8j [24[ [25j [26]. 
This analysis requires a more careful study and is work in progress. 
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